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Abstract 



i-^ I We develop a theory describing the effects of many-particle Coulomb cor- 

["t I ■ relations on the coherent ultrafast nonlinear optical response of semiconduc- 

f^ . tors and metals. Our approach is based on a mapping of the nonlinear op- 

tical response of the "bare" system onto the linear response of a "dressed" 
^ , system. The latter is characterized by effective time-dependent optical tran- 

^D ' sition matrix elements, electron/hole dispersions, and interaction potentials, 

^s^ ■ which in undoped semiconductors are determined by the single-exciton and 

CN ■ two-exciton Green functions in the absence of optical fields. This mapping 

^^ . is achieved by eliminating the optically-induced charge fluctuations from the 

Q^ \ Hamiltonian using a Van Vleck canonical transformation. It takes into ac- 

count all many-body contributions up to a given order in the optical fields as 
well as important Coulomb-induced quantum dynamics to all orders in the 
optical field. Our approach allows us to distinguish between optical nonlin- 
earities of different origins and provides a physically-intuitive interpretation 

Q \ of their manifestations in ultrafast coherent nonlinear optical spectroscopy. 

O 
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I. INTRODUCTION 

In semiconductors and metals, the Coulomb interactions among electrons and holes 
strongly affect the optical spectra close to the onset of absorption.lil In undoped semi- 
conductors, the e-h attraction leads to exciton resonances, while in doped quantum wells 
and metals it leads to the Fermi-Edge- Singular ity.H In the vicinity of the resonances, the 
interactions between photoexcited e-h pairs and Fermi sea electrons also affect strongly 
the time-dependence of the e-h polarization measured in ultrafast nonlinear optics exper- 
iments. This allows one to obtain new information on many-body correlations, not easily 
accessible with other experimental techniques, by using transient pump-probe and wave- 
mixing spectroscopy,^ The interpretation of such experiments, however, requires solving a 
challenging time-dependent many-body problem. Indeed, during femtosecond time scales, 
the Coulomb interactions dominate the time-dependence of the experimentally-measurable 
polarization. In the coherent limit, the standard quasi-equilibrium approaches no longer 
apply. This led to the development of alternative approximations for undoped semiconduc- 
tors. The most important of those is the Semiconductor Bloch Equations (SBE's) approach, 
which treats exciton-exciton interaction effects on the nonlinear optical response within the 
Hartree-Fock approximation.aO 

Recent experiments, however, revealed pronounced features in the amplitude and phase 
dynamics of the nonlinear polarizations that could not be captured by the SBE's.Q In Ref. 
7, for example, pronounced features in the phase and amplitude dynamics of the four- 
wave-mixing polarization were attributed to non-Markovian dephasing due to the Coulomb 
correlations, neglected within the Hartree-Fock approximation. In Ref. 8, it was shown that 
the Boltzmann picture of relaxation is inadequate for extremely short optical pulse durations. 
In Ref. 9, it was observed that the time decay and spectral width of the four-wave mixing 
polarization did not correlate with the decay of the time-integrated signal, as would have 
been expected based on the SBE's. Experimental evidence for the effects of the bound 
biexciton state (not described by the SBE's) was presented, e.g., in Ref. 10. Correlation 
effects due to the repulsive exciton-exciton interaction were found to significantly enhance 
the four-wave-mixing signal for negative time delays as function of magnetic field, even in 
the absence of bound biexciton states.B Excitation-induced dephasing effects were observed, 
e.g., in Ref. 12. 

Experiments such as those mentioned above point out the need to incorporate in the the- 
ory exciton-exciton correlations beyond the time-dependent Hartree-Fock approximation. 
This problem has been addressed within the equations-of-motion approach for density ma- 
trices, which, to a given order in the optical fields, leads to a closed system of differential 
equations that determine the nonlinear polarizations.BEJ'EJ The solution of such equations, 
however, must address the formidable four-body problem of two electrons interacting with 
two holes. With the exception of one-dimensional systemsEj, the determination of such biex- 
citon states is numerically very demanding. Various approximation schemes, devised mainly 
for calculations of the third-order polarization, x^^\ have been recently proposed.l3't3'EJ The 
e-h Coulomb correlations also strongly affect the coherent polarization femtosecond dynam- 
ics in doped quantum wells.ll^O Such effects originate from the dynamical Fermi sea response 
to the photoexcited valence holes,p'0 and can be observed for frequencies close to the Fermi 
enegy, where the absorption spectrum is dominated by the Fermi-edge singularity and the 



e-e quasiparticle scattering is inhibited.O'ta'EJ The ultrafast dynamics of the Fermi-edge 
singularity was addressed in Ref. 16. 

In this paper, we develop a new approach for describing correlation effects in ultrafast 
nonlinear spectroscopy. Our approach is motivated by the fact that the nonlinear polar- 
ization measured in pump-probe or four-wave-mixing experiments arises from the linear 
response to the probe optical field (propagating along direction k^) of a system described by 
the Hamiltonian H + Hp{t), where H is the two-band interacting Hamiltonian in the absence 
of optical fields and Hp{t) describes the coupling to the pump optical field (propagating along 
direction kp). Terms nonlinear in the probe optical field do not contribute to the third-order 
polarizations measured along directions k,- (pump-probe) or 2kp — k^ (wave-mixing). They 
can also be neglected when the probe is weaker than the pump field, as is often the case ex- 
perimentally. However, unlike the "bare" Hamiltonians H, the Hamiltonian H + Hp{t) does 
not conserve the number of e-h pairs. The idea is to "block-diagonalize" the Hamiltonian 
H + Hp{t) in the basis of states with fixed number of e-h excitations, such as, e.g., exciton 
states, biexciton states etc. This can be achieved to any given order in the optical fields 
by using Van Vleck canonical transformations,Eifc2l generalized to the time-dependent case. 
This way we obtain an effective time-dependent Hamiltonian, Heff{t) which, does conserve 
the number of e-h excitations. Note that for the time-independent case, such a procedure 
is analogous to the Schrieffer- Wolff transformation used to eliminate the charge fluctuations 
in the Anderson Hamiltonian,c3o where Hqs corresponds to the Kondo Hamiltonian. 

The above procedure allows us to cast the experimentally-measurable nonlinear polar- 
izations in the form of the linear polarization of an effective two-band system. The latter 
is characterized by time-dependent band dispersion relations, interaction potentials, and 
transition matrix elements, which can be obtained to any given order in the pump optical 
field. In particular, as we will show, to second order in the pump field the above effective 
parameters can be expressed in terms of the one- and two-exciton Green functions for the 
"bare" Hamiltonian H. Such a formulation allows us to interpret the various dynamical 
features in the nonlinear absorption spectra, originating from the correlation effects, within 
the familiar framework developed for linear spectroscopy. 

In particular, in linear absorption, the strength and spectral position of the excitonic 
resonances are determined by (i) the dipole transition matrix elements and (ii) the time 
evolution of the e-h pair wavefunction. Here we establish a similar picture for the nonlinear 
absorption spectrum, which applies for any pulse duration: (i) The e-h pair photoexcited 
by the probe interacts with the e-h pairs photoexcited by the pump via the Pauli principle 
(Phase Space Filling effect), static exciton-exciton interactions (Hartree-Fock contribution), 
and exciton-exciton scattering processes that involve the exchange of center-of-mass mo- 
mentum (correlation effects). Such effects are described by a renormalized transition matrix 
element, (ii) Similar interactions also lead to a pump-induced renormalization of the semi- 
conductor band energies and e-h interaction potential, which affect the time evolution of 
the (probe-induced) e-h pair wavefunction. The corresponding time evolution operator that 
creates the excitonic resonance is governed by the Hermitean effective Hamiltonian H^ff(t) 
which, importantly, has the same form as H. Therefore, the nonlinear absorption spec- 
trum can be viewed as arising from the probe-induced optical transitions to renormalized 
("dressed") exciton states, which are determined by the Schrodinger equation with Hamil- 
tonian Hes{t). One, thus, should distinguish between two physically-distinct origins of the 



optical nonlinearity: optically-induced renorinalizations of the transition matrix elements 
and time-dependent parameter renormalizations in the Hamiltonian. The latter not only 
shift the bound state energy, as would have been the case within the two-level system ap- 
proximation, but also changes its wavefunction and Bohr radius that determine the exciton 
oscillator strength. Importantly, due to the Coulomb interactions, Hes{t) does not commute 
with itself at different times (unlike for monochromatic photoexcitation), which leads to 
quantum dynamics and memory effects that, as we show, are not captured by the third- 
order nonlinear response ( x^^"*) or by using few-level system truncations. Here we provide 
a systematic method for including all the above effects in a physically-intuitive fashion. 

The paper is organized as follows. In Section II we set up the problem and briefly outline 
our final results. In Section III we describe the canonical trasformations that eliminate all 
optically-induced charge fluctuations up to the second order in the pump field. In Section 
IV we obtain the linear response of the system to the probe optical field. In Section V 
we focus on the pump-probe nonlinear polarization and identify the effective Hamiltonian 
and optical transition operator. In Section VI we derive the four-wave-mixing polarization. 
Section VII concludes the paper. 

II. PROBLEM SET UP 

In the absence of optical fields, the excitation spectrum of a two-band semiconductor is 
described by the Hamiltonian 

^ = E <^-q ^-q + E« + Eg) a^ a^ + Ke + ^M + Kh, (l) 

q q 

where a], is the creation operator of a conduction electron with energy e'i and mass me, 
6Lq is the creation operator of a valence hole with energy e^ and mass rrih, Vee, Veh, and 
Vhh describe the e-e, e-h, and h-h interactions, respectively, and Eg is the bandgap. For 
simplicity, we consider a single parabolic valence band and suppress the spin degrees of 
freedom. 

Ultrafast pump-probe and wave-mixing spectroscopy measures the response of such a 
system to photoexcitation by a pump optical field, £'p(t)e*^p'""*'^''* (with Ep{t) centered 
at t = 0), and a probe optical field, Er(t)e'^^'''~^'^p^^~'^'' (with Erit) centered at t = r), 
propagating along directions kp and k^, respectivelv. and have central frequency Up and 
time-delay r. In the rotating wave approximation,B3 the coupling to the optical field is 
described by the Hamiltonian Hp(t) (pump), 

Hp{t) = dpit) [e^'^p-'-f/t + H.c] , (2) 

and the Hamiltonian H^-lt) (probe), 

Hrit) = dr{t) ^e^^^-^+i'-priji + H_c.] , (3) 

where W = SqCtq^-q is the "bare" optical transition operator, di{t) = fiEi{t),i = p,T are 
the dipole energies, and /i is the interband transition matrix element. The total Hamiltonian, 
including the coupling to the optical fields, is 



Htotit) =H + Hpit) + Hrit) - hujp J2 «li «q- (4) 

In the following, we will absorb the last term into H. 

The pump-probe and four-wave-mixing signals are determined by the polarization 

P(t)=^e-^-p*(vl/(t)|?7|M/(t)), (5) 

where the state |^E'(t)) satisfies the time-dependent Schrodinger equation with the Hamil- 
tonian Htotit). In the following, we distinguish between pump and probe optical fields by 
concentrating on the polarizations that propagate along direction k^ (pump-probe) and 
2kp — k^ (four-wave-mixing). For a weak probe, the nonlinear polarizations arise from the 
linear response to the probe-induced perturbation Hr{t) of the coupled-two-band system 
system described by the time-dependent Hamiltonian H + Hp{t). Note that, for x^^^ along 
the above dierctions, the above is also true even if the probe and pump fields have com- 
parable strength. In contrast to H, this Hamiltonian does not conserve the number of e-h 
pairs, so we transform it into an effective Hamiltonian Hes{t) that has this property by using 
canonical (Van-Vleck) transformations.ElrE3 

Specifically, the Hilbert space of the bare semiconductor states (in the absence of opti- 
cal fields) consists of disconnected subspaces, which can be labelled by the number of e-h 
airs. In this basis, H has block-diagonal form. The Hamiltonian Hp[t) couples the differ- 
ent subspaces by creating/ annihilating e-h pairs. A unitary transformation that "block- 
diagonalizes" the Hamiltonian H + Hp{t), e.g., up to the second order in pump field, has 
the form e~^^e~^^[H + Hp{t)]e'^'^e^^, where the antihermitian operators Si{t) and S2{t) cre- 
ate/annihilate one and two e-h pairs, respectively. In order to achieve this in higher orders 
in the pump field, we should also include in Si(t) and 5*2 (t) processes that create/annihilate 
odd and even numbers of e-h pairs, respectively. Correspondingly, the state |\I^(t)) is given 
by 

|^(t)) = e^iWe^2W|<l>(t)) (6) 

where all the effects of the probe optical field are contained in the state |$(t)). After 
the transformation operators Sl"' (t) are calculated and the effective Hamiltonian Hcsit) 
is derived, the polarization Eq. (^ with |\I^(t)) given by Eq. (||) can be found in the first 
order in the probe field by identifying all contributions that propagate along directions k,- 
(pump-probe) or 2kp — hr (four-wave-mixing). 

Let us summarize our results for undoped semiconductors. The pump-probe nonlinear 
polarization, P^^ can be presented as a sum of two terms, Pk^{t) = P^^'^ (t) + P,^^''^'"^ (t), 

where P^^'^ (t) and Pj^ ^'^^'^'(t) are polarizations due to the transitions from the full valence 
band ground state to exciton and biexciton states, respectively. In particular, 

PiT\t) = -^e^*^— ^(*-^)/_*^rft'rf.(t')(0|t/c.c(t)W(t,t')f/eU(t')|0), (7) 

where 

UUtm=T.M,.{t)aibl^\0) (8) 



is the effective optical transition operator to exciton states, Mk(t) being the correspond- 
ing matrix element. Here U{t,t') is the time-evolution operator satisfying the Schrodinger 
equation 

th^V({t,t') = H,s{t)U{t,t'), (9) 

where Hee{t) is the effective semiconductor Hamiltonian that acts on the single e-h pair 
subspace. We show that Hcg{t) can be cast in the same form as H: 

i/effW = E<W«i«q + E<W^-q^-q- E Veh{k,k';t)albl^b.^, a^> . (10) 

q q kk'q 

Explicit expressions for the renormalized matrix elements and parameters of Hef[{t) in terms 
of the single- and two-exciton correlation functions are derived in Section V. The biexcitonic 
contribution to the pump-probe nonlinear polarization is given by an expression similar to 
Eq. ([7p, with the corresponding two-pair transition operator Ul^^^^{t) (see Section V). 

Note that in the absence of pump-induced photoexcitations, Ep{t) = 0, the familiar form 
of the linear absorption spectrumcJ is simply recovered from Eq. (^ by setting U{t,t') — > 
e^i^(*^* ) and Ul^^{t) -^ f/^. The similar second-quantized forms of Hamiltonians Hcs{t) and 
H suggest an interpretation of the excitonic effects in the nonlinear absorption spectrum 
in terms of probe-induced optical transitions to "dressed" exciton states. This analogy 
becomes more clear if one expresses the time-evolution operator U{t, t') in the basis of states 

i^{t,t')=j:\^^\t)){¥:\t% (11) 

n 

which satisfy the Schrodinger equation for the Hermitean Hamiltonian Hes{t) with the initial 
condition that, before the pump pulse, |$^^^(t)) coincide with the exciton eigenstates of H. 
As we will show in the next Section, the time-dependence of the parameters in the 
effective Hamiltonian Eq. ([lOD follows that of the pump amplitude. Therefore, as can be 
seen from Eq. (|^), one can replace U(t, t') by e~*^(*~* )/^ when the probe comes at time t' = t 
after the pump. Thus the optically-induced corrections in Hcs{t) affect the pump-probe 
spectra only for negative probe delays or when the two pulses overlap in time. Furthermore, 
their effects are most pronounced in pump-probe spectroscopy, where, unlike in four-wave- 
mixing, they already contribute to the third-order nonlinear polarization. On the other 
hand, the optically-induced renormalization in the matrix element Mk(t), which includes 
contributions from the Phase Space Filling as well as from the exciton-exciton interactions 
(see Section V), is finite even after the pump optical field is gone and therefore dominates 
the pump-probe signal for positive probe delays. The same is true for the biexcitonic 
contribution to the pump-probe nonlinear polarization, which, as discussed below, vanishes 
for negative time delays. 



III. PHOTOEXCITATION BY THE PUMP OPTICAL FIELD 



A. Elimination of single pair transitions 



Our goal is to "block-diagonalize" the Hamiltonian H+Hp{t) by performing the canonical 
transformation Eq. (^. In this subsection we determine the transformation operator Si{t) 
that eliminates the single-pair pump-induced transitions. We start with a slightly more 
general problem, useful for calculating higher order nonlinear polarizations, by deriving a 
differential equation for the operator S{t) which ensures that all odd-pair charge fluctuations 
are eliminated from the transformed Hamiltonian. Ei 

In the time-dependent Schrodinger equation [ihd/dt — Htotit)]]"^ (t)) = we substitute 
|\l/(t)) = e'^''*^|x(t)) and, after acting with the operator e""^^*-* on its Ihs, we get 



-Sit) 



dt 



;5(*)|^(t)) = e-^W [Hp{t) + H,{t)] e^W|x(t)). 



(12) 



Since both S{t) and Hp{t) create/annihilate odd numbers of e-h pairs and H conserves the 
number of e-h pairs, all odd-pair contributions to |x(^)) originate from terms in the Taylor 
expansion of the exponentials in Eq. (|12[) of the order of (5*) "'*' (i^^ — Hj or (S) ^ Hp, n 
being an integer. Following Ref. 22, we determine S{t) by requiring that all such terms, 
which describe pump-induced charge fluctuations, cancel out. 

Using the Baker-Campbell-Hausdorff expansion for an arbitrary operator A{t), 



-s(t)A(t)e'(^) 



coshS{t) + smhS{t) A{t) 



(13) 



where S{t) is a superoperator defined by S{t)A(t) = [A{t), S{t)], we find that S{t) should 
satisfy 



ihsmh[S{t)]-^ = smh[S{t)]H + cosh[S {t)]Hp{t). 



(14) 



Such a form is inconvenient, however, since the Ihs includes multiple commutators of S{t) 
with its time-derivatives. These can be eliminated by applying S{t) sinh~^ S{t) (defined by 
its Taylor expansion) to both sides of Eq. (|1^. Using the relation S-^ = ^, we obtain 



^^d^ = [H, S{t)] + S{t) coth[S{t)]Hp{t), 



(15) 



with the initial condition S{—oo) = 0. 

Eliminating all contributions to Eq. (^) with odd numbers of e-h pairs by using Eq. 
(p!5|), we obtain after collecting the remaining terms 



cosh[5(t)] 



d 



\X 



it)) = smh[Sit)] Hp{t) \x{t)) + e'^Wi7.(t)e^W|x(t)). (16) 



This equation, however, is not particularly useful because it does not have the form of 
a Schrodinger equation. After some manipulations involving adding and subtracting the 

operator ih^^ — H and inserting the identity iSj 5* = 1, it can be rewritten in the form 



;4-^)w^)) 



'cosh[S{t)] - l] [S{t)Y' S{t) U^^ - h\ \x{t)) 
+ smh[S{t)]H,{tMt)) + e-^«i/.(t)e^W |x(t)). 



(17) 



We now simplify the first term on the rhs of Eq. (|T^) by making use of Eq. (13). After some 
algebra, we finally obtain 



ih 



dixit)) 
dt 



H + tanh 



S{t) 



H,{t) 



|x(t)) + e-^Wif.(t)e^W|x(t)). 



(18) 



Eq. (^) now has the form of the Schrodinger equation for |x(t)), with a time-dependent 
Hamiltonian depending on S{t). The above derivation generalizes the results of Ref. 22 to 
time-dependent systems. 

The operator 5*1, which eliminates the single-pair excitations, can be obtained by lineariz- 
ing Eq. ([I5|). The nonlinear in 5*1 terms in Eqs. ([T5|) and ( plSf ) describe the corrections in the 
effective semiconductor parameters of the order of E^ or higher. Using the decomposition 
of the anti-hermitian operator S'i(t), 



Si(t) = Si(t)e- 



T\{t)e 



(19) 



where S|(t) and Si(t) create and annihilate single e-h pairs, respectively, we obtain the 
following equation for S|(t): 



ih 



asl(t) 



dt 



/f,SI(t)l -d^{t)U\ 



which has the formal solution 

Sl(t) ='-[' rft'rfp(t')e*^(*'-*)/^[/te-*^(*'-*)/^ 

it J —CO 



(20) 



(21) 



Note that, since the Hamiltonian H conserves the number of e-h pairs and the optical tran- 
sition operator U^ creates a single e-h pair, S|(t) also creates a single e-h pair. Furthermore, 
since both H and U^ conserve momentum, so does S|(t). 

It is convenient to represent S|(t) as a sum of one-body and many-body parts. 



where 



sl(t) = st(t) + sL(t), 



T.i{t) = Y.PM.t)<bl, 



(22) 



(23) 



describes the photoexcitation of a single e-h pair (or exciton) with zero total momentum, 
while the many-body operator S^^(t) contains quartic and higher order terms in the single- 
particle operators. As we discuss later, in undoped semiconductors, S^(t) describes the 
exciton-exciton interactions. After substituting Eq. (^) into Eq. (pO|) and collecting all 
terms quadratic in the single particle operators, we obtain 



ih^ P.(q, t) = [m + 4 + e:;) P.(q, t) - ^ i;(q - q')Px(q', t) - dp{t), 



(24) 



where v{q) is the Coulomb potential and AQ = Eg — hup is the detuning. Thus, the 
amplitude Px{q^,t) is just the linear polarization due to the pump field. This can also be 
seen by substituting Eq. (|^) into the polarization equation Eq. (|^) (with S^^ = 0) in the 
absence of the probe field. 

Note that, due to the interaction terms in H, the commutator in the rhs of Eq. (pOD 
generates the following terms, quartic in the single-particle operators: 



ULit) = E ^(q) [^x(k, t) - P.(k - q, t)] UUk, q) 

kq 



with 



ULik, q) = E (4-q^lk4'+qak' 



«k-q^-k^-k'+q^-k' 



(25) 



(26) 



Thus, we get the following equation for T,]^^{t): 



ih 



d^Ut) 



dt 



H, SL(i)l + UUt) 



which has the formal solution 



with 



^L(i) = ^ E ^(q) f dt' [P.(k - q, t') - P.(k, t')] Ulik, q, t - t') 



f/t,(k,q,t) = e-^^*/^f/t,(k,q)e^^*/^ 



(27) 



(28) 



(29) 



As can be seen from Eq. (^), S]^^(t) describes the photoexcitation of a single e-h pair 
accompanied by Coulomb interactions with additional photoexcited or Fermi sea carriers. 
Note that the time evolution of S^^(t) is driven by the linear pump-induced polarization 
and thus, in contrast to Sx(t), only indirectly depends on the pump amplitude dp(t). 

We now turn to the Schrodinger equation Eq. (0). The leading-order corrections to 



the effective Hamiltonian are given by the terms linear in 5"!. With the help of Eq. (19) we 
obtain 



d 

ih— Hcsit) 

dt ^ ' 



|^(t)) = _^ fe^^'^p'- 



(e2^'^pr[f/t,Sl(t)]+/7.c.)|x(t)) 
+e-^iW/f,(t)e^iW|x(t)). 



where 



H^f^it) = H - ^ {[uMit) 



+ H.C. 



(30) 



(31) 



is a time-dependent effective Hamiltonian that conserves the number of e-/i pairs and Si(t) 
is given by Eq. (^Oj) with initial condition S][(— oo) = 0. The first term in the rhs of 



Eq. (|30|) describes the pump-induced two-pair creation/annihilation processes. These can 
be ehminated by performing a second canonical tranformation, as described in the next 
subsection. 

Concluding this subsection, let us address the condition of validity of our approach. It 
is useful to write down a formal solution of Eq. ( pOD in the basis of the N-hole many-body 
eigenstates, \aN) with energies EaN, of the Hamiltonian H. Here, for example, A^=0 gives 
the semiconductor ground state |0), |al) denotes the exciton eigenstates, \a2) denotes the 
biexciton eigenstates, etc. In this basis, the solution of Eq. pO|) can be written as 

{PN + l\J:\{t)\aN) 



{pN + l\W\aN) 



Tl J —oo 



where we separated out the detuning Afi and denoted Ai?^^ = -E/jat+i — EaN] the width F 
describes the effects of dephasing processes not included in H (e.g., due to phonons). It can 
be seen that for resonant excitations (small Afi) the rhs is of the order of dptp/h. Thus, for 
short pulses, this parameter justifies the expansion in terms of optical fields. Note that, for 
off-resonant excitation, this expansion is valid even for longer pulse durations. 

B. Elimination of t^vo pair transitions 

We now perform the second canonical transformation in order to eliminate the first term 
in the rhs of Eq. (^), which describes biexcitonic transitions. We define 

|x(t)) = e^^W|<l>(t)), (33) 

where, in the absence of the probe field, \^{t)) has fixed number of e-h pairs. Again, using 
the anti-Hermicity of 5*2 (t), we decompose it as 

S2{t) = S2(t)e-2*''p-'" - J:l{t)e^'^^-\ (34) 

where S2(t) and S2(t) create and annihilate two e-h pairs, respectively. Substituting Eq. 
(p3|) into Eq. (|30| ) and requiring that all two-pair terms cancel out, we obtain similar to the 
previous subsection the following equation for Tj2{t) to the leading order, 

^^^=[HMit)]-i^[i:\it),u% (35) 

which has the solution 

St(t) = ±- t dt%(t')e*^(*'-*)/"[Sl(t),f/t]e-^^(*'-*)/^ (36) 

As we will show later, T^Kt) includes all exciton-exciton interaction contributions to the 
third-order four-wave-mixing polarization but only affects the pump-probe signal via higher 
order {E^) parameter renormalizations. 
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IV. PROBE INDUCED OPTICAL TRANSITIONS 

In the previous section, we "block-diagonalized" the Hamiltonian H + Hp{t) (up to the 
second order in the pump field) in the basis of states which, in the absence of the probe 
optical field, conserve the number of e-h pairs. In the presence of the probe field, the state 
|$(t)) (Eq. (|33D) satisfies the Schrodinger equation 



ih^^\^{t)) = H^simit)) + dr{t) [e^^-'-^-^^^^fZ/W + H.c] |<l>(t)), (37) 

where Hf,s(t) is given by Eq. (|3TD and 

[/t(t) = e-^^We-^^W[/te^^We^^W (38) 

is the (transformed) optical transition operator. Following the standard procedure, we obtain 
in the first order in the probe field 

|$(t)) = |i _ 1 /■* d,{t')U{t,t') p''-'-+*-''^f/;(t') + H.c] I |<l>(°)(i')), (39) 

where U{t, t') is the time-evolution operator satisfying 

tn^U{t,t')=H,s{t)U{t,t'), (40) 

and |$*^°)(t)) = U{t, —oo)\0) is the time-evolved ground state |0). Since Hes{t) conserves 
the number of e-h pairs, |$*^*'^(t)) contains no e-h pairs. In undoped semiconductors, the 
only such state is the ground state so that \^^^\t)) oc |0). 

Using Eq. (|3|) for |$(t)), the polarization P(t) = /ie-*^^*(^(t)|f/|^(t)) = 
/xe-^'^''*($(t)|f/r(t)|$(t)) can be written as 

P{t) = --^e-'^^' f dr{t') 

tl J — OO 

-(0| [e**^- '■+'"''"[// (t') + e-*'^-'-'"''"f/r(t')] Uit',t)U,it)\0)}. (41) 

Explicit expressions for U^{t) are derived in the next sections. 



V. PUMP-PROBE SPECTRUM 

A. Pump probe polarization 

In this section we derive the closed form expression for the nonlinear polarization prop- 
agating in the direction of the probe optical field. In order to extract the pump-probe 
polarization from Eq. (|4l|) , one should retain only contributions that are proportional to 
e^^^'^ and are independent of the pump direction kp. To the leading order, the transition 
operator U^ (t) can be obtained by expanding Eq. (|38|) in terms of ^i and S2 (given by Eqs. 



(|T^) and (|34D ) using the Baker-Campbell-Hausdorff formula, Eq. 
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UKt) = UUt) + e-^^^^U^oit) + f/liWe-^^^''^ + Ul,,{t)e^'^^^ + Ul{t)e'^^^\ 



(42) 



(Note that the higher-order corrections do not contribute also to x measured along the 
probe direction.) Here the operator 



uLit) = u^ + ^[^iit),[u\^\{t) 



+ 



Sl(t),[f/t,Si(t) 



creates a single e-h pair; the operator 

t/o^(t)=[f/t,Si(t) 
conserves the number of e-h pairs; the operator 

creates two e-h pairs; the operator 



uUt) = \u\^2{t) 



+ 2 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



while all other terms annihilate the ground state |0). Substituting Eqs. ( ^Sf ) and (^Of ) into Eq. 
(PI) and retaining only terms proportional to e^^^^, we obtain -Pk^(t) = Pj^'j^^'^ (t) + Pj^^'^™ (t), 
where 



annihilates a single e-h pair; and the operator 

1 



f/t,Si(t)l,Si(t) 



f/t,sl(t)l,sl(t) 



ul{t) = [j:i{t),u[ 

creates three e-h pairs. Furthermore, it is easy to see that 

Ur{t)\0) = [e'^^'^Uoit] + U.iit)e^'^''''] 
and 

Uhtm = \uUt) + e-^'-^^U^.it) + Ul^Me'''^' + uUt)e'^'''''] |0), 






and 



-J (biexc) 
k. 



it) 



H^ ilir-r-iujpit-T) 

' h 



rft'rf.(t')(0|f/e.,(t)W(t,t')f/L(t')|0), (50) 



dt'd^{t'){0\Uue.cm{t,t')Ul,,,{t')\0). (51) 



In the derivation, all terms containing Uq cancelled each other out. Note that the above 
expressions apply to both undoped and doped semiconductors. 

We have thus expressed the nonlinear polarizations, Eqs. (^)-(|5T|), in terms of the linear 
response to the probe optical field of a system described by a Hermitean time-dependent 
effective Hamiltonian, Eq. (|3T|), that does not include any charge fluctutations. The first 
term, Eq. (^OD, descibes the excitonic (single-hole state) contribution to the pump-probe 
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polarization. Eq. (|50|) contributes to the pump-probe spectrum at frequencies close to the 
exciton energy. Since ^biexc(^) creates two e-h pairs, the second term, Eq. (^l]), describes 
a biexcitonic contribution, not included in the SBE's. Note that S|(t), and hence Ul^^^^{t) 
vanishes if the probe pulse arrives before the pump pulse. Therefore P^^ '*^™ (t) is unim- 
portant for negative time delays, for which the time-dependence of Hes(t) contributes the 
most. Biexciton (two-hole) states also affect P^^'^ (t) by acting as intermediate states in 
the renormalization of the parameters entering into the optical transition operator Uexdt) 
and the effective Hamiltonian ifeff(^); as described later. 

In the following sections, we focus on undoped semiconductors. We show that f/exc(^) 
and Hf,s{t) have the same operator form as their "bare" counterparts U and H, respectively; 
they only differ in the parameters. This fact allows us to view the non-linear polarization 
as arising from the linear response to a probe-induced optical transition from the ground 
state to a nonstationary "dressed" exciton state evolving in time with Hamiltonian Hcs{t). 

B. Optical transition matrix elements 

In this subsection we derive explicit expressions for the effective transition matrix ele- 
ments. We start by showing that, when acting on the ground state |0), as in Eq. (|5D|), the 
operator U^x^{t) has the same operator form as W. Substituting E{(t) = Sx(t) + S^^(t) 
into Eq. ( ^31) and using the second-quantized expression for ^^(t), Eq. (|23|), we obtain after 
some algebra 

uLMo) = E [i - l^x(k,t)|2] aibl^io) - 1 [Sx(t)SL(t) -i/.c] f/t|0) 

-iSxx(t)SL(t)f/^|0). (52) 

In deriving Eq. (^) we used the fact that Si(t) annihilates the ground state |0) and that 

SL(t)|0) = Sxx(t)f/t|0) = Sxx(t)St(t)|0) = (53) 



The latter relations can be obtained from Eqs. ( PSJ ) and ( PB| ) by noting that H does not 
change the number of e-h pairs. 

It can be seen from Eqs. (|5^) that, since both S^xl^) ^^^ ^i(^) create a single e-h pair 
[see Eqs. (^), (0), and (|28|)], the transition operator U^-^^it) also creates a single e-h pair. 
Since, in addition, the operators Sj^(t) and S]^(t) conserve the momentum, as discussed 
above, we have 

f/]xcW|0)=E^k(t)4&lk|0). (54) 

k 

By comparing to the "bare" transition operator, W\0) = Z^k^^k^lklO)) ^^ ^^^ ^^^^ ^^^ 
parameter 

Mk(t) = (0|6_kakt/eL(t)|0), (55) 

should be interpreted as an effective time-dependent optical transition matrix element. 
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To derive an explicit expression for Mk(t), it is convenient to work in the exciton basis. 
We introduce the operators Bj^ and -Bn,q that create/annihilate an exciton with center-of- 
mass momentum q and energy i?„q in the relative motion state 0„(q): 

Bl^ = j:Mk + (3c0aU^bl^, (3= 1^ (56) 

with rrie and rrih being electron and hole mass, respectively, and define the corresponding 
exciton contribution to the linear polarization as 

^.w = E^x(k,t)0:(k). (57) 

k 

The effects of the exciton-exciton interactions can be conveniently described by introducing 
the correlation function 

S..(nWq'; n; t) = ^{0\B.^,^,B^^,^^>^Ut)Bi,\0), (58) 

which represents the contribution to the probabihty amplitude of pump-induced transi- 
tions between single-exciton, 5„q|0), and two-exciton, i?^,_q,i?^,q, |0), states due to the 
exciton-exciton interactions. Such transitions also occur in the absence of interactions, with 
amplitude given by Eq. (^) but with S]^(t) replaced by S|^(t). This determines the first 
(Phase Space Filling) term in the rhs of Eq. ([5^). As we discuss below, by restricting in 
Eq. (|58|) to two-exciton states with zero center-of-mass momentum q' = 0, one recovers the 
static exciton-exciton interaction contribution included in the SBE's. 

Substituting alb_^^ = J2n4'nO^)Bnfi into Eqs. ( p4D and ( p2D and inserting the closure 
relation into the last term of Eq. (|52|) we obtain from Eq. (|5^) that Mk(t) = M^ (t) + 
M^\t), where M^^^\t) = 1 - |Px(k,t)|^ is the usual Phase Space Filling and M^''\t) 
comes from the exciton-exciton interactions: 

nn'm' 

-$:^;'W0n'(k)Sxx(riWO;n;t)] 

~Q Yl 0m(k)$*S*x(?^Wq';m;t)Sxx(^Wq';n;t), (59) 

O III 

where $„ = I]k0n(k). The first term in Eq. (^) contains both static exciton-exciton in- 
teractions (Hartree-Fock contribution) and exciton-exciton scattering processes (correlation 
contribution). It describes a process where the photoexcitation of an e-h pair by the probe 
optical field is accompanied by the pump-induced creation and subsequent annihilation of 
an e-/i pair with zero center-of-mass momentum. Correlations contribute to this term via 
the process during which the pump photoexcites an e-h pair with finite momentum due to 
scattering with the probe-induced e-h pair, which then scatters back to zero momentum 
and subsequently annihilates. The last term in Eq. (|59|) describes exciton-exciton scatter- 
ing processes (not included in the SBE's) where the pump photoexcites an e-h pair with 
finite momentum due to scattering with the probe-induced e-h pair, which subsequently 
annihilates by giving its momentum back to the probe-induced e-h pair. 
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The contribution of the biexciton transitions, |0) -^ ^biexc(''')|0)' ^^ ^^^ pump-probe 
polarization is given by Eq. (^). The biexciton optical transition operator f/biexc(^)' given 
by Eq. (|^), can be obtained in a similar manner by inserting the closure relation expressed 
in the basis of momentum eigenstates and noting that only two-pair states contribute. The 
result reads 

UL.citm = l E K^Un'm'c^;n;t)Bl^,Bl_^,\0), (60) 

n'm'q'n 

where Tjyr^{n'm'q';n;t) is given by Eq. (pH]). Below we show that T,^n'm'q';n;t) can be 
expressed in terms of the two-exciton Green function for the bare Hamiltonian H. 



C. Effective Hamiltonian 

We now derive the second-quantized form of the effective Hamiltonian Hcs{t), whose 
formal expression is given by Eq. (^Tj). Here we primarily focus on the excitonic contribution 
to the nonlinear polarization and therefore, according to Eq. (|50|) , we only need to consider 
how Hesit) acts on single-pair states. One has to distinguish between two contributions to 
the pump-induced correction to the Hamiltonian, | f [f/, S|] + H.c.j, which come from the 

decomposition Si(t) = S^(t) + Sj^^^(t). The first contribution, due to the one-body single- 
exciton operator S]^(t), is also present in two-level systems. Using Eq. (|23|) , we obtain after 
straightforward algebra 

^ {[U, St] + H.c.) = E ReP.(q', t) - ^ ReP.(q', t) {al,a^, + bl^,b^^,) . (61) 

q' q' 

The second term in the rhs represents a correction to the electron and hole band energies 
(the first term plays no role and can be dropped). This is the origin of the energy shifts 
(ac-Stark effect) well-known from the dressed atom picture.B3'El Here, however, the shifts 
are momentum-dependent. This leads, in particular, to an additional time-dependent cor- 
rection to the exciton Bohr radius and wavefunction absent within the two-level system 
approximation.ll3 

To derive the many-body correction to Hcs{t), described by St^(t), we note that, ac- 
cording to Eqs. (1^) and (^4|), the time evolution operator U(t,t'), and hence Hesit), only 
act on the state U^^^,{t)\0) that describes a single e-h pair with zero momentum. Therefore, 
it is sufficient to describe the action of Hes{t) on the basis of states alj)_^\0). Since both 
St^(t) and f/t conserve momentum and create a single e-/ipair, we obtain that, when acting 
on states alb^_k\0), 

UJ:Ut) + H.c. = E f^eh(p, k; t)a;6lp6_^ak, (62) 

kp 



where 



f/e/,(p,k,t) = (0|6_pap \UJ:Ut) + H.c.] 4&-k|0). (63) 



Following the steps described in the previous subsection and transforming to the exciton 
basis, Ueh{p,^,t) can be presented as 
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[/e;.(p, k, t) = Y. 0n'(p)C(k) [$™'S..(nWO; n; t) + $:;,,S:,(nm'0; n'; t)] , (64) 

nn'rra' 

with Hx^ln'm'O] n; t) given by Eq. (|58|) . Putting everything together, we finally arrive at the 
following second-quantized expression for Hc^it): 

Hes{t) = E<(^)«l«q + E<W^-q^-q - E t^eh(k,k';t)4ak' &lk&-k', (65) 

q q kfc'q 

where 

<W = 4 + c?p(t)ReP,(q, t), t = c,v, (66) 

are the time- dependent renormalized dispersion relations of electrons and holes, respectively, 
and 

v,h{k, k'; t)=v{k- k') - dp{t)U,h{k, k'; t) (67) 

is the (nonlocal) effective e-h interaction. 

It should be emphasized that the pump-induced parameter renormalizations in Hesit) 
and U^^^it) affect the polarization Eq. (|50|) in different ways. Indeed, the dependence of 
Hesit) on time follows that of the pump pulse: the last terms in Eqs. ( |66|) and (|67|) are 
proportional to dp{t). Therefore, they do not affect the time-evolution operator U{t,t') 
when the probe comes at times t' after the pump pulse, but do contribute to the pump- 
probe signal for negative probe delays or when the probe overlaps in time with the pump 
pulse. On the other hand, the transition operator Eq. (^31) depends on the pump only via 
S|(t) and hence, according to Eq. (|21|), is determined by the pump amplitude at all past 
times t' < t. For positive time delays, the dynamics of the pump-probe spectra is therefore 
governed by the renormalized transition matrix element and the biexcitonic contribution Eq. 



(pT|). In this case, the latter is determined by the time-evolution operator corresponding to 
the bare Hamiltonian H. 

We now discuss the role of the pump-induced corrections in the effective Hamiltonian 
Hes(t). As discussed in the Introduction, the spectral position and strength of the exciton 
resonance is determined by the time evolution operator U{t, t'). In the linear response theory, 
this operator satisfies Eq. (^Of ) with Hamiltonian H, while in nonlinear absorption the latter 
must be replaced by the effective Hamiltonian Hcs{t). Therefore, Hes{t) plays an important 
role in determining the properties of the excitonic resonance observed in the pump-probe 
spectrum. First, the time-dependent change in the band dispersion relations in Hes{t) leads 
to a shift in the position of the excitonic resonance, which is, of course, the origin of the 
ac-Stark effect. Second, in the presence of interactions, Hcs{t) does not commute with itself 
at different times, which leads to memory effects. Importantly, such effects only arise to 
second or higher order in the Hamiltonian correction dp{t)AH{t) = —^dp{t) ([U, S|] -|- H.c.) 
and therefore do not contribute to the third-order polarizations. To illustrate this point, let 
us expand U{t, t') in terms of the optical field using the Magnus expansion:E2l 



W(t,t')=e-*^(*-*')/"exp 
where 
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.n=l ^'■ 



(68) 



r .H(t-,-t') .Hiu-t') 

Ui{t,t') = I dtidp{ti)e'^^ AH{ti)e~'^^ , (69) 



U2{t,t')= f dp{ti)dti r dt2dp{t2 



.H(u-t'} .H(t-,-t') .HHy-t') .H(ty-t') 



(70) 



etc. In the absence of Coulomb interactions, [H, AH{t)] = and only the first term survives. 
Ui (t, t') leads to the energy shifts if one does not expand the exponential. However, in order 
to describe the Coulomb-induced quantum dynamics emerging from the second and higher 
order terms in Eq. (|68D, one cannot restrict the pump-probe polarization to the third order 
in the optical fields, but rather must use the full U(t,t'). Below we describe a simple 
physically-intuitive way, based on the similarity between Hesit) and H, for calculating the 
response to the probe in all orders in the optically-induced terms in Hcs{t). 

We start with the pump-probe polarization Eq. (|^). At time t', the probe photoexcites 
the system into a single-pair state U^^^{t')\0) , Eq. (15^) , with the momentum wavefunction 
Mk(t'). Note that, for negative time delays, Si(t') = and hence Mk(t') = 1. The subse- 
quent evolution of such state up to time t is governed by Hes(t) and describes the formation 
of the exciton resonance. The momentum wavefunction of the time-evolved state, $k(^), 
satisfies the Schrodinger equation 

^^9$k(t^ = Kit) + slit) + rfp(t)?7eh(k, k, t)] $k(i; t') - J2 ^ehiK k', t)$k'(i; t'), (71) 

with the initial condition $k(^';^') = Afk(^')- Correspondingly, the polarization P^^^'^ (t), 
Eq. (pUj), takes the form 

ptr\t) = f dt'd^{t')Y^Ar^m^{t-t') (72) 

Thus, polarization dephasing arises from two different pump-induced effects: the change 
in the transition matrix element Mk(t) and the time-evolution of the state $k(^;^') due to 
the time-dependent coefficients in the Schrodinger equation Eq. ([7I|). Note here that, for 
Ep{t) = 0, Mk(t) = 1 and Eq. (|72D reduces to the familiar Fermi's Golden Rule expression. 
Note also that Mk(t) and the effective Hamiltonian parameters are time-independent in the 
case of monochromatic photoexcitation. 

The Hamiltonian Hes{t) thus determines new ("dressed") exciton states, which depend 
on the time-dependent optical field amplitude and, as in the linear response theory, can be 
obtained by solving the Schrodinger equation Eq. ([7T|) . For example, this point of view pro- 
vides a natural explanation of why, for below-resonant photoexcitation, the pump-induced 
coupling to higher exciton states leads to a rigid exciton blue-shift instead of bleaching, as 
would have been the case if only the Is bound state were included. l^'oa'^H To see this, note 
that, to first approximation, the renormalization of the band dispersion relations, Eq. (|66|), 
results in an increase of the e-h effective mass. This in turn leads to an optically-induced de- 
crease in the Bohr radius a^, which, as is known from linear absorption, directly determines 
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the strength of the exciton resonance. Since the exciton strength is inversely proportional 
to a% {d is the dimension), a decrease in Bohr radius enhances the exciton strength and 
competes with the Phase Space Filling reduction. Our approach also allows us to study the 
effect of the pump on the exciton strength for short pulse durations, where the adiabatic 
limit does not apply, and compare it to the Fermi Edge Singularity. t3'cllll3 This and other 
applications will be discussed elsewhere. 

D. Exciton exciton interactions 

In this section, we obtain the equation for the correlation function 'L^^{n'm!al\ n; t) which 
determines the renormalization of the parameters in the effective Hamiltonian and the tran- 
sition matrix elements. 

As can be seen from Eq. (|58D , the time dependence of Ti^n'm'q'; n; t) is determined by 
the operator S]^^(t). In order to understand which processes contribute to T,xx{'n'm'q_';n;t), 
it is useful to present T,]^^{t), defined by Eq. (pSD, in second-quantized form: 

"""'" ° „. J,„,„, («+l)!«!(rn+l)!rn!^""«- P^' 

where the operator 

SL(t)= E S„„[{p';k';k;p};t]aJ,,---aJ,, 6Lk;---&Lk' , 
{pp'kk'} 

x6-k™---&-kiap„---api (74) 

describes the photoexcitation of a single e-h pair accompanied by its Coulomb interactions 
with n additional photoexcited or Fermi sea electrons and m additional photoexcited valence 
holes. 

Substituting Eqs. (|7^ and (^) into Eq. (^), one can see that Sjj^(t) couples only to 
^n-imit) s-iid T,l^^_i(t) but not to terms with higher n or m. This defines an iterative scheme 
for obtaining successively higher amplitudes. As a result, in undoped semiconductors, only 
a finite number of low-order terms S]^^(t) contribute to the nonlinear optical response to a 
given order in the optical fields. For example, only the following interaction terms contribute 
to the third-order nonlinear polarizations: S|Q(t) (describing interactions between the pho- 
toexcited e-h pair and an additional electron), ^li{t) (describing interactions between the 
photoexcited e-h pair and an additional hole), and S|]^(t) (describing interactions between 
the photoexcited e-h pair and an additional electron and hole). Thus, a calculation of the 
third-order nonlinear polarization requires the solution of a four-body problem. Higher- 
order interaction terms do not contribute because the corresponding Sjj„(t) give zero when 
acting on states with up to two e-h pairs. 

We now derive the equation for J]^x{n'm'q'; n; t), which is more convenient for describing 
the above processes in the case of excitons. Starting from the definition Eq. (|58D , using Eq. 
(p8|) for Sj^x(t), and going to the basis of exciton eigenstates of H, we obtain after some 
algebra 



ll J —CO 

X Y. ^(m'n'q'; m"n"q"; t - t')t^(q") F„„. (q") 

n"m"q" 

X 5: [P.(k', t') - P.(k' - q", t')] ct>l„{\^ - /3q") (75) 

k' 

where 

Km(q) = E i</'n(k + q) - 0„(k)] 0;;,(k + /?q) (76) 

k 

and G{m'n'c(] m"n"c^'; t) is the two-exciton Green functional defined as 

Gim'n'q!; m"n"q"; t) = {0\Br,^^,Bm',-^'Bl„_^„{t)Bl,^^„{t)\0), (77) 

with i?^ q(t) = e~*-'^*i?^ qC*^*/^. Thus, the problem of calculating the effective parameters 
is reduced to the two-exciton Green function Eq. ([77|) with the "bare" Hamiltonian H. The 
Dyson equation for G can be obtained in a straightforward way starting from the equations 
of motion for the exciton Heisenberg operators 

G(mVq'; mnq; t)= e-^(^".-<^+^'"-)*(0|5„,,q,fi„,,_q,<_q5t q|0) 

ih Jo 
X E K"Z"(q"-q)G(mVq';mV'q";t'), (78) 

m,",n",q" 

where Kl'm'(q) = 'i^q-Pnn'(q)-^mm'(— q) IS the exciton-exciton interaction potential. 

The Hartree-Fock approximation corresponds to retaining only the first term in the rhs of 
Eq. (1781). Substituting the latter into Eq. (^) and using our previous results for the effective 
transition matrix element and Hamiltonian parameters, one obtains the contribution to the 
nonlinear polarization that is of first order in the exciton-exciton interaction. It can be seen 
that, in this case, only two-exciton states with zero center-of-mass momentum contribute 
to the nonlinear polarizations. The correlation effects due to exciton-exciton scattering 
involving the exchange of center-of -mass momentum are described by the second term 
in the rhs of Eq. (|75|). To lowest order in the exciton-exciton interaction, these can be 
obtained by solving Eq. ([78|) iteratively. This is equivalent to a second- or higher-order 
Born approximation to the self-energy that determines the effects of the exciton-exciton 
interactions on the nonlinear optical polarization. In the case where opposite-spin excitons 
are photoexcited, Eq. (|78D can also be used to study the effects of bound biexciton states. 

As can be seen from Eq. ([75|), the two-exciton Green function acts as a memory kernel 
that determines the effective parameters at time t in terms of the linear polarization at 
earlier times. If one only retains the first term on the rhs of Eq. ( |78D (Hartree-Fock approx- 
imation), such memory effects decay exponentially as determined by the dephasing time. 
The significance of the correlation term in Eq. ([78| ) is that, under certain conditions, the 
memory kernel can become non-exponential (non-Markovian memory effects). The impor- 
tance of such effects for short pulse excitation will be discussed elsewhere. Note here that in 
the case of monochromatic excitation, for which S^ is time-independent, the correlations 
simply change the magnitude of the various parameters. 
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VI. FOUR WAVE MIXING POLARIZATION 

Let us now turn to the polarization along the four- wave-mixing direction, ks = 2kp — k,-. 
The latter can be obtained by substituting Eqs. ( ^S]) and (^) into Eq. ( pD and retaining 
only terms proportional to e*''^'". After some algebra, we obtain 

PFWM(t) = !|^e*k3r-.,(*-.) f' dr{t'){0\[UU{t',t)Ul^y,{t) - {t ^ t')]\0), (79) 

where 

Ulwuit) = UWlit) - St (t)t/St (t), (80) 



and 



wl,{t) = \[j:l{t)f-j:l{t), (81) 

is the "irreducible" biexciton operator satisfying 

^^^^PpWlO) = HWlitm - ^f/tst(t)|0). (82) 

In the above, we used the fact the St^(t)|0) = 0. Thus, the four wave mixing polarization 
is determined by Sx(t), which describes the linear pump-induced polarization, and the 
biexciton transition operator Wpp{t), which includes all exciton-exciton interaction effects 
in the four-wave-mixing signal. One can express VFjJ (t)|0) in terms of the two-exciton 
Green function Eq. (^) by following the procedure of the previous section. Note here that, 
unlike for the pump — probe polarization, the pump-induced corrections to the Hamiltonian 
do not contribute to the third-order polarization along the direction ks. Therefore, the x^^^ 
calculations of the four-wave-mixing signal neglect the optically-induced change in the time 
evolution operator U{t,t'), while similar pump-probe calculations include this effect only to 
first order in perturbation theory. As discussed in the previous section, such calculations 
do not properly treat the shift of the exciton resonance and neglect the memory effects 
and quantum dynamics induced by the interactions in Hes{t). The latter are especially 
important in the case of the Fermi-edge singularity.!!^ 

Finally, let us discuss the connection between our approach and other works. In Refs. 13 
and 14, it was shown that the third-order nonlinear polarization can be expressed exclusively 
in terms of single-pair and two-pair creation operators, corresponding to our St(t) and 
W^ (t). Similarly, one can obtain x^^^ along the four- wave-mixing direction by simply 
replacing H^sit) in Eq. ([79|) by H. The connection to Ref. 14 can be established by noting 
that (a2|14^jjp(t)|0) coincides with the lowest term in the expansion of the Hubbard operator 
|0)(a2|. Furthermore, by projecting Eq. (p2D into two-pair states, one can identify the 
amplitude (0|6 ii.6_k"3k'-qak+qW^pp(i^)|0) with the four-point density matrix introduced by 
Axt and Stahl.li^l This connection is less obvious for our pump-probe results, Eqs. ( |50D 
and (0). The reason is that, as discussed above, these include important contributions 
beyond x^^^- However, by only keeping the lowest term in a perturbative expansion of the 
time-evolution operator [i.e., by restricting to x^^A^ one can obtain after straightforward but 
laborious algebra the expression for the pump-probe polarization exclusively in terms of 
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E|(t)|0) = M* dt'd,{t')e-''"^'~''^/^U^\0), (83) 

/i J — OO 

and 

W,\{t)\0) = 1 f dt'diity-'"'^'-''^''' f/ts|(t')|0), (84) 

ri J —CO 

with i,j=p, r. The corresponding expression is rather tedious and is omitted here. By 
comparing Eq. (|H3|) to Eqs. (^ and (|2|), it can be seen that St(t)|0) = St(t)|0). It can 
also be seen from Eq. (|8^) that Wjjp(t)|0) coincides with the solution of Eq. ([821) . The third- 
order polarization is thus determined by the states pyjj(t)|0), which, for i ^ j, describe the 
effects of the interference between the pump and the probe. 



VII. CONCLUSIONS 

In conclusion, we discussed a general theoretical approach for describing Coulomb many- 
body effects in the ultrafast coherent dynamics in semiconductors and metals. Our approach 
provides a complete description of the correlation and non-Markovian memory effects by 
reducing the calculation of the nonlinear optical response of the "bare" system onto the linear 
response of a "dressed" system described by the effective time-dependent Hamiltonian, Eq. 
(|65|). Furthermore jt allows us to study the coherent ultrafast nonlinear optical response of 
Fermi sea systems.li3 

As is known,El the calculations of the third-order nonlinear polarizations starting from 
equations of motion must address the time evolution of exciton and biexciton states. In 
view of the unresolved numerical problems facing a straightforward calculation of biexcitonic 
effects (except in one dimension), our approach has the advantage of providing a physically- 
intuitive picture for viewing their manifestations in the ultrafast nonlinear optical response. 
Close to the exciton frequency, these effects can be taken into account via the renormalization 
on the electron/hole dispersions and interaction potential as well as the probe-induced 
optical transition matrix elements. Our approach also allows us to treat the Coulomb- 
induced quantum dynamics arising from the semiconductor response to such perturbations, 
which cannot be captured within the x*^^^ approximation. Note that, for monochromatic 
excitation, the exciton-exciton correlations simply change the magnitude of the effective 
parameters, while for ultrashort pulses they also affect their time-dependence over times 
comparable to the dephasing time. 

Our theory separates the effects due to the renormalization of the electron/hole dis- 
persions and interaction potential from the renormalization of the probe-induced optical 
transition matrix elements. These affect the nonlinear polarization in a different way, as can 
already be seen in the absence of Coulomb interactions by recalling the well-known "dressed 
atom" results for a two-level system.c3 In the latter case, the shift in the energies results in an 
overall resonance shift (e.g., ac-Stark effect) which cannot be correctly described within x^^^ 
calculations. Such effects are described here by the time evolution operator for the effective 
Hamiltonian ifeff(^)- In semiconductors, the additional quantum dynamics (not captured 
by X^^^) arising from the non-commutativity of the effective Hamiltonian at different times 
can be simply described with the Schrodinger equation for the effective exciton Hamiltonian 
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Hes{t). On the other hand, the renormahzation of the transition matrix elements, due to 
Phase Space Filhng and, in semiconductors, exciton-exciton interaction effects, leads to an 
overall factor in the polarization that can be calculated to third-order in the optical fields. 
Thus our theory identifies a number of competing physical processes, each leading to 
different nonlinear dynamics and spectral features. In undoped semiconductors, we can ob- 
tain the parameters determining the strength of the above nonlinearities in terms of the 
"bare" semiconductor one- and two-exciton Green functions. Finally, we also established a 
connection between our approach and x^^^ calculations in undoped semiconductors. We be- 
lieve that the present theory is well-suited for interpreting the physics conveyed by ultrafast 
coherent nonlinear optics experiments in semiconductors and metals using pulses as short 
as 10 fs. 
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